arXiv:astro-ph/0512195vl 7 Dec 2005 


Scale-invariant Perturbations from Chaotic Inflation 


Christian T. Byrnes and David Wands 

Institute of Cosmology and Gravitation, University of Portsmouth , Portsmouth POl 2EG, United Kingdom 

(Dated: December 7, 2005) 

Vacuum fluctuations in the inflaton field driving chaotic inflation with a quadratic potential give 
a red spectrum of primordial density perturbations, n = 0.97. However angular fluctuations in 
an 0(N)-symmetric quadratic potential have a very nearly scale-invariant spectrum, n = 0.9998. 

We investigate the possibility that these isocurvature field perturbations could give the dominant 
contribution to the primordial density perturbation after inflation. 
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I. INTRODUCTION 

Probably the simplest model for inflation in the very early universe is an accelerated expansion driven by an inflaton 
field with a quadratic potential [lj. In this paper we consider the spectrum of primordial perturbations in the general 
case of a quadratic potential 


V = ^m 2 \<j> 2 \, ( 1 ) 

where (f> is an IV-component scalar field with rotational symmetry. In this case the radial field \<f>\ plays the role of 
the inflaton field and orthogonal fluctuations are isocurvature perturbations during inflation as they do not effect the 
energy-momentum tensor to linear order. For example, a massive real scalar field corresponds to N = 1, whereas a 
massive complex scalar field corresponds to the case N = 2, and the potential energy is independent of the phase of 
the field. The model 0 appears as a special case of inflation driven by multiple massive fields, previously studied in 
models of double inflation [2, Q and assisted inflation 0. 

Vacuum fluctuations in the inflaton field yield a lower bound on the amplitude of the primordial density perturbation. 
This spectrum deviates from exact scale-invariance, n = 1, due to the slope of the potential and effective mass of 
the inflaton field, characterised by the dimensionless slow-roll parameters e and r) a respectively. In general slow- 
roll inflation the scale-dependence is given by n — 1 ~ — 6e + 2ri a {jjjj. For a quadratic inflaton potential we have 
e = rj a ~ 0.01 and the perturbations have a red spectrum with n ~ 0.97. Combined with the tensor-scalar ratio, 
r ~ 16e ~ 0.13, this is marginally compatible with current observations Q. 

But a primordial density perturbation can also be generated from vacuum fluctuations in fields other than the 
inflaton. Examples are the curvaton proposal Q, modulated reheating <§[ ||] or asymmetric preheating in a nearly 
symmetric potential E3- The general form for the scale-dependence of perturbations in an isocurvature field \ 
during slow-roll inflation is given by — 2e + 2ij x 0, 12]. Thus for inflation driven by the quadratic potential 0 with 
e = r] a = rj x the isocurvature perturbations have a very nearly scale-invariant spectrum. 

This paper is organised as follows. In the next section we briefly review the dynamics of the homogeneous inflaton 
field driving chaotic inflation, and then calculate the scale-dependence of the radial and angular field perturbations. 
We find the spectral tilt of the isocurvature perturbations is second-order in the slow-roll approximation. In Section 
III we then investigate a possible mechanism that might transfer the angular fluctuations to the primordial density 
perturbation and compare this contribution with that from inflaton perturbations. We conclude in section IV. 


II. INFLATION 


We consider a scalar model with effective action 


S = 



M 2 1 

E1r- 


( 2 ) 


where Gjj is the scalar field metric. We will assume the potential, V, has a minimum at 4>k = 0 and that as it 
approaches the minimum it can be described by a symmetric quadratic potential of the form given in Eq. 0 where 


l</> 2 | = ■ 

I,J 


( 3 ) 
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Henceforth we take Gij = 6ij so that the fields are canonically normalised real fields. 

We can exploit the rotational symmetry to set y = </>j = 0 initially for all J > 2 and take a = <j>\ to be the inflaton 
field. Thus we can write 

V = \m 2 a 2 , (4) 

along the classical trajectory for any number of real scalar fields in a symmetric quadratic potential. This inflaton 
field obeys the standard equation of motion for a single scalar field 

a + 3 Ha + m 2 a = 0 . (5) 

X = 0 remains the classical solution for the spatially homogeneous background field. However there will be quantum 
vacuum fluctuations in both the radial field, Sc r, and orthogonal fluctuations, <5y. Radial fluctuations in the inflaton 
field, <7, yield adiabatic density perturbations on large scales, while any fluctuations away from y = 0 correspond 
to non-adiabatic fluctuations. Non-adiabatic perturbations can in general affect the density perturbation on large 
scales, but our symmetric potential has a vanishing first derivative orthogonal to the inflaton, V x = 0, and hence the 
fluctuations Sx are strictly isocurvature perturbations and do not couple to the adiabatic perturbations at first-order 

13. 


A. Slow-roll parameters 


We conventionally define a dimensionless slow-roll parameter 


1 


f= Mh(VA* 

167T \ V ) ’ 


( 6 ) 


where V„ = dV/da 1 which describes the slope of the potential relative to the Planck scale. Two more slow-roll 
parameters describe the curvature of the potential with respect to the inflaton field (down the potential) and with 
respect to the orthogonal isocurvature field perturbations 


Va 


Vx 


M 2 , V aa 
8tt V ’ 
^pi Vxx 
8t r V ' 


(7) 

( 8 ) 


For a general potential there are additional parameters at lowest order in the slow-roll approximation proportional to 
V ax [T]J, but these are zero due to our assumption of rotational symmetry. 

At next order in the slow-roll expansion we also define 


2 _ Afp; 

Scr — 




2 _ Afp; V a V aXX 


(8tt) 2 V 2 

For the symmetric quadratic potential 0 we have 

e = Va=Vx = 


r = 

Sy — 


(8tt) 2 V 2 


47TfJ 2 


^ 2 = £ 2 = o . 


(9) 


( 10 ) 


Note that the slow-roll parameters are thus independent of the mass, m, and depend solely on the value of the inflaton 
field relative to the Planck scale. 

During slow-roll inflation driven by a quadratic potential the integrated expansion, or number of e-folds from the 
end of inflation, is given by 


N = 

Hence the slow-roll parameter m is given by 


H dt ~ 


27r(a 2 - a 2 ) 
Mp 2 


1 ( 27T(7g \ 

£ ~ 2N\ + NM 2 J ’ 


( 11 ) 


( 12 ) 


where e e = Mp^/Ana 2 ~ 1 at the end of inflation. Comoving scales comparable to our present horizon size crossed 
Hubble scale during inflation when N cm b — 60 Hi and hence we have e cm b — 0.008 at that time. 
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B. Inflaton perturbations 


Perturbations in the inflaton field are gauge-dependent, but are most easily described in the spatially flat gauge 
jlfij . Inflaton perturbations with comoving wavenumber k on spatially flat slices obey the wave equation [l6j 


5a + 3H5a + 


■V m 


87t d 
MpjO 3 dt 


a 3 a 2 

H 


5a = 0. 


This can be succinctly written in terms of the variable u = a5a E3 


k 2 -) u = 0, 

z 


(13) 


(14) 


where a prime denotes derivatives with respect to conformal time dr = dt/a and z = aa/H. To leading order in the 
slow-roll approximation we have z" jz ~ (2 + 9e — 3 t)<j)/t 2 . 

In order to calculate the power spectrum to leading order in the slow-roll approximation it is sufficient to match 
the solution for vacuum fluctuations on small scales, k/a^> H , 


1 


u = 


—ikr 


\/2k 


and match it to the growing mode solution on large scales, k/a <C H, 

u (X z , 


(15) 


(16) 


(which corresponds to a constant curvature perturbation H5a/a ) at the time when that mode crosses the Hubble 
scale: k = aH. We thus obtain the power spectrum for inflaton field perturbations on large scales 


_ 47rfc 3 9| 



H H' 

\h) 

a 27r 


(17) 


k—aH 


The scale-dependence of the spectrum ra follows from the time-dependence of the quantities in the square-bracket 
evaluated when k = aH. The spectral tilt is then given by 


_ dlnVa r , 0 
An CT = ~ —6e + 2 r] a . 

ainfc 

Hence for a quadratic potential with e = ri a ~ 0.008 we have 

A n a ~ —0.03 . 

Note that the running of the inflaton spectral index is given by 

d&n a _ 2 


din k 


- -24e" + 16 er, a - , 


(18) 

(19) 

( 20 ) 


where £ 2 = 0 for our quadratic potential and hence we have a a — — 8e 2 and hence of higher-order in the slow-roll 
parameter. 


C. Isocurvature perturbations 


Field fluctuations, 5\, orthogonal to the inflaton field are automatically gauge-independent [l3j since y = 0 in 
the homogeneous background solution. These are isocurvature and hence they obey the wave equation for field 
perturbations in an unperturbed FRW metric: 


5 X + 3H5 X + l^ + V xx )5 X = 0. 


a* 

In conformal time this can be written in terms of the variable v = a5x as 

„// > 


v 


a 


aV xx -1 v = 0 


a 


( 21 ) 


( 22 ) 
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where to leading order in slow-roll a"/a ~ (2 + 3e)/r 2 and m 2 a 2 ~ 3e/r 2 . Hence the equation above can be rewritten 
as 

v = 0 where v = ^ + 0(e 2 ). (23) 

Note that for isocurvature perturbations in a quadratic potential the terms O(e) in Eq. 12211 cancel out. Thus the 
equation of motion |23|) to first-order in slow-roll is the same as that for a massless scalar field in de Sitter (although 
the actual field is neither massless, nor in de Sitter), v is a constant up to second-order in slow roll parameters and 
hence the approximate solution on timescales of order the Hubble time, is given by Hankel functions of order v — 3/2. 
Matching the general solution to the small-scale asymptotic vacuum gives 



—ikr 


\/2k 



(24) 


To obtain a more accurate solution over many Hubble times (without using the slow-roll approximation), on 
large scales, we note that the equation of motion for isocurvature perturbations is exactly the same as that for the 
homogeneous inflaton field m in the long-wavelength limit where gradient terms can be neglected. Therefore we have 
Sx = A(k)a on large scales. Matching this solution with the growing mode of Eq. (I2H) at Hubble exit, i.e. at k = aH 
or equivalently at kr = — (1 + e) leads to 


V x {k) 


47rfc 3 
( 2^)3 


I Yf * 


(1 — e)H 


1 2 


2na 


k—aH 


We can use the first-order slow-roll solution 


at Hubble-exit, to obtain 



47T 

T 



V x {k) 



( m \ 2 

[l - —el 

37T 

\Alpi) 

3 J 


(25) 


(26) 


(27) 


<r 2 /M|> ; ~ N/2tt gives the overall time-dependence of the spectrum during inflation. The only scale-dependent term 
in this equation is the slow-roll parameter, e, evaluated at Hubble-exit. 

It is straightforward to calculate the spectral tilt 


A n x = 


d In V x 
din k 



(28) 


Since this is second order in slow roll it is very small, A n x = —0.0002. The running of the isocurvature spectral index 
is even smaller, being third-order in slow roll 


dAn x 40 c. 
= din A; ~ ~T £ ' 


(29) 


In Appendix El we show that our calculation of the spectrum and its scale dependence can be confirmed using a 
more accurate Hankel function solution around Hubble exit, including the second order slow-roll correction to v in 

Eq. 01. 


III. PRIMORDIAL DENSITY PERTURBATIONS 

We have shown that inflation driven by a quadratic potential produces an almost scale-invariant spectrum of 
field perturbations on super-Hubble scales by the end of inflation. To give rise to the observed anisotropies in the 
microwave background sky and large-scale structure in our Universe today, these field perturbations during inflation 
must produce density perturbations in the radiation dominated era after inflation. These primordial perturbations 
are usefully characterised in terms of the dimensionless density perturbation on spatially-flat hypersurfaces. For linear 
perturbations we define 

HSp 


P 


(30) 
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This is equivalent to the perturbed expansion, 5N = HSt , up to a uniform-density hypersurface some time after 
inflation. 

Because we are interested in perturbations on scales very much larger than the Hubble scale, the local expansion, 
N, is given in terms of the background solution for the local values of the scalar fields on an initial spatially-flat 
hypersurface during inflation. This is known as the “separate universes” approach E3- Thus for one or more scalar 
fields, </>/, during inflation, we have |1S| 


(3D 

Lyth and Rodriguez j2jj have recently pointed out that the extension of this result to second-order also allows one 
to calculate the non-Gaussianity of primordial perturbations due to non-linear dependence of the expansion after 
Hubble exit on the initial field values. Non-Gaussianity from adiabatic field fluctuations in single-field inflation are 
small (first-order in slow-roll parameters |2ll l22jl but provide an important constraint on non-adiabatic fluctuations 
in multi-field models (2^, [3 |2f|. 

The power spectrum of primordial density perturbations can thus be written as a sum of contributions from the 
power spectra of individual fields 



(32) 


This can then be split into a contribution from inflaton perturbations, 5c r, and the orthogonal isocurvature perturba¬ 
tions during inflation 





(33) 


The existence of a very nearly scale-invariant spectrum of isocurvature perturbations on large scales during inflation 
thus allows the possibility of producing a very nearly scale-invariant spectrum of primordial density perturbations 
after inflation. These perturbations are on wavelengths far greater than the Hubble scale at the end of inflation and 
so their subsequent evolution should be very well described by the long-wavelength approximation in which spatial 
gradients can be neglected and hence is independent of wavenumber. Thus any conversion of isocurvature to curvature 
perturbations after inflation will be independent of scale and yield a spectral tilt n — 1 = A n x . 


A. Perturbations during inflation 


Field perturbations in the inflaton field during inflation correspond to perturbations along the background trajectory 
and thus the perturbed expansion can be calculated anytime after Hubble-exit, to give 


Pc.inf = 



Pa- 


Using Eq. in and the slow-roll approximation GS> we have 


Pc,inf 


167T 

TOCT 2 

2 

4 

( nn N 

) 2 

[ Y | 27RT el 

3 

YK J 

k=aH 

37T 

^ Mpi , 

) 

[ + m;U 


(34) 


(35) 


Given the WMAP normalisation, P^- ~ 2 x 10 -9 [26}, then for m ~ 10 -6 Mpi fluctuations in the inflaton field produce 
primordial density perturbations of the observed magnitude. On the other hand for m < 10 -6 Mpi the inflaton 
perturbations are below the observational bound and we would require an additional contribution from isocurvature 
perturbations to generate the primordial density perturbations. 

We also see that primordial density perturbations due to inflaton fluctuations will have the same spectral tilt <P|) 
as the inflaton perturbations on large scales during inflation, conventionally written as n = 1 + An ff ~ 0.97. 

It has been shown that the non-Gaussianity from inflaton fluctuations during slow-roll will be negligible with 
non-linearity parameter f^L — — 5e/3 (2lj . 
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B. Perturbations from instant preheating 

Perturbations in the orthogonal fields do not affect the expansion history during inflation due to our assumption 
that the potential, V , is symmetric under rotations. But these isocurvature perturbations during inflation can give 
an additional contribution to the density perturbation in the radiation era if the symmetry of the self-interaction 
potential CD is broken by other interactions of the field. In particular they will produce density perturbations if the 
local expansion, N, is sensitive to the initial isocurvature perturbations, dN/d\ 7 ^ 0 in Eq. (1331) . 

We will consider a model of instant preheating at the end of inflation where the inflaton field’s energy density can 
be transferred to a preheat field during the first oscillation of the inflaton field [27j • This simplifies the calculation of 
particle production as a function of the initial field values, making it possible to estimate the effect on preheating of 
the isocurvature perturbations during inflation. It may also be the case that instant preheating maximises the effect 
of the isocurvature perturbations which might be dissipated if reheating took place over many oscillations. 

In this scenario the preheat particles i/j are created through an interaction term g 2 cfr ip 2 / 2 , during a brief period 
around <j> = 0, when their effective mass vanishes. The energy density in the y-field is then “fattened” by the coupling 
to the inflaton field as <f> rolls back up the potential because of their effective mass, m ^ = g 2 4> 2 . Through a Yukawa 
interaction A ipuiUj the preheat field can then decay into to particles, decaying most rapidly when the ^ particles reach 
their maximum effective mass, as the inflaton field reaches its maximum. Depending on the coupling constants this 
process may be so efficient that all further decay of the inflaton field can be neglected [ 23 , US ■ 

The efficiency of instant preheating is sensitive to the phase of a complex scalar field if the real and imaginary 
components (<pi and ^> 2 ) have different couplings to the preheat field. We consider an interaction term 

Ant = ~{gl4>\ + g&l)^ 2 , (36) 

which explicitly breaks the rotational symmetry of the inflationary potential CD- 

Particle creation in the preheat field first occurs when the adiabatic condition fails, \rh^\ = m , shortly before 
the inflaton first passes through its minimum. Denoting this time with a so that e.g. cr* = er(t*), the adiabatic 
condition is broken when 


I o'* I = | ct*| 2 g 


(37) 


where the effective coupling constant is given by 

g 2 {9) = g\ cos 2 9 + g\ sin 2 9 (38) 

and <f >i = cr cos 9 and <j >2 = <rsin0. In practice |<r*| ~ mMpi and is independent of g, and hence independent of the 
phase 9, for m <C gMp\. 

The time interval during which particle creation takes place is 

At* - = i< 7 *r l/a r 1/2 - (39) 

After the inflaton field has passed through the origin for the first time the occupation number of the tj.’ field with 
wavenumber k is j27| 


rife = exp[—7r(fcAf*) 2 ] . 

We integrate rife to give the total number density of y-particles produced 

1 

n v = —cr 

x (2 tt) 3 


d 3 kn.k = (27rAt*) 


-3 


(40) 


(41) 


The ratio of the energy density of the y-field to that of the inflaton field at the point when the inflaton would be 
at its maximum after the first explosive creation of i/j particles is p^/p a ~ g 5 / 2 (10~ 6 Mp\/m )i / 2 [ 2 ^ Thus if g > 1 
we see that most of the inflaton’s energy density is transferred to the preheat field during the first oscillation for 
m ~ 10 -6 Mpi. This constraint on g is dependent on to and the energy transfer is more efficient for to <C 10 -6 Mpi, 
which would correspond to the inflaton perturbation being smaller than the observed primordial perturbation. 

Assuming the decay products are massive, non-relativistic particles, then the expansion from an initial spatially- 
flat hypersurface to a final uniform-density hypersurface is given by N oc (l/3)logn^,. Hence the linear metric 
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perturbation I® is given by £ = (1/3Thus if the effective coupling, g, is dependent upon the phase of the 
complex field we have 


where 


_ 1 d\ogg _ F(.R,tan6>) 

~ 2 86 ~ 2 cr 


F(R, tan 6) 


tan 0(1 — R 2 ) 
R 2 + tan 2 6 


and R=— . 

92 


(42) 


(43) 


We only need to consider the case R < 1 and 0 < 8 < 7r/2, due to the symmetries of this model. The power spectrum 
is then given from Eqs (E3) and CS by 


F 2 (R, tan 9) / m \ 2 
~ 12tt \m7i) 


(44) 


Comparing this with Eq. 13511 we see that either V^ t i so or V^.inf can dominate the primordial density perturbation 
depending on the values of R and 9. The asymmetric preheating creates the dominant primordial perturbation if 


F(R, tan 9) > AN ~ 240. 


(45) 


We require R < 1/480 for the isocurvature perturbations to dominate in any range of 9. If R <C 1/480 then the 
isocurvature perturbations dominate for 240I? 2 < 9 < 1/240, while larger R will have a smaller range of suitable 9. 
The range of 9 where isocurvature perturbations dominate will be larger if IV < 60. We require that the inflaton 
trajectory is almost exactly along fa, the field with the smaller coupling constant g\. Although such small values for 
9 may, a priori, seem unreasonably small, it may be that the interaction cmt plays a role in selecting the inflaton 
trajectory during inflation, but such an investigation of the initial conditions for the inflaton are beyond the scope of 
this paper. 

The WMAP normalisation that — 2 x 10 -9 constrains the three free parameters m, R and 9. Equation 1351) 
shows that the primordial perturbations from the inflaton alone are of the required size if m ~ 10 _6 Mpp But 
if m 10 ~ e Mpi then we require some additional source of primordial density perturbations, and perturbations 
produced at preheating from isocurvature field fluctuations are of the required size if 

F(R, tan 8) = 200 ( 10 6m pi \ (46) 

\ m J 

Using the formalism of |2C| and identifying the non-linear expansion as N ex (1/3) log {n^(8)), and assu ming the 
isocurvature field perturbations themselves are Gaussian, the non-Gaussianity of the primordial perturbation (221 ] can 
be given as 


_ 5N, gg _ 5 (R 2 — tan 2 9) 

tNL ~ 6 N 2 ~ 3(1 -R 2 ) sin 2 9 

This is of order unity, and hence likely to be unmeasurable, for 9 ~ R, but can become large for smaller 9. For example 
if R = 10~ 3 , then the perturbations from asymmetric preheating dominate in the range 0.00026 < 9 < 0.0039 and 
f jvl lias its maximum deviation from zero of —24 at 9 = 0.00026, but is far smaller for most of the allowed range of 9. 


IV. DISCUSSION 

We have shown that isocurvature field perturbations during chaotic inflation with a symmetric, quadratic potential 
have an almost exactly scale-invariant spectrum. We have shown that the deviation from scale-invariance is second- 
order in slow-roll parameters. The rotational symmetry of the inflationary potential leads to the scale-invariance of 
the isocurvature spectrum and keeps these field perturbations decoupled from the adiabatic density perturbations in 
the inflaton field driving inflation. But the spectrum of initial isocurvature perturbations can be transferred to the 
primordial density perturbation if interactions after inflation break the rotational symmetry. The spectral index of 
the generated primordial density perturbation is then given by n ~ 0.9998. 

The relative importance of adiabatic inflaton fluctuations versus isocurvature field fluctuations in determining the 
primordial density perturbation in the radiation-dominated era long after inflation depends upon the perturbation in 










the local expansion |nj , represented by the dimensionless perturbation £ in Eq. <0. In contrast with the inflaton, the 
amplitude of the density perturbations resulting from orthogonal perturbations is dependent upon the physics after 
inflation. Given that the inflaton and isocurvature field perturbations on the Hubble scale near the end of inflation 
are of similar magnitude, and that the inflaton spectrum is red so grows on larger scales, we actually require 


ON 

dx 


> 



(48) 


for the primordial perturbations due to fluctuations in an orthogonal field to dominate over the conventional inflaton 
result. 

It turns out to be rather difficult to construct a simple model that can produce primordial density perturbations 
from the isocurvature perturbations that are larger than those from adiabatic perturbations in the inflaton field in our 
chaotic inflation model with a symmetric potential. The expansion during slow-roll inflation is sensitive to the value 
of the inflaton field, dN/da oc e -1 / 2 and hence the primordial power spectrum is boosted by a factor e -1 relative 
to the inflaton field perturbations. If primordial density perturbations due to isocurvature field perturbations are to 
dominate, then their effect on the expansion history must be boosted by a larger factor. 

In the simplest case of a single complex field, with a U(l) symmetric potential, one could simply postulate that the 
efficiency of reheating at the end of inflation is an arbitrary function of the phase, 9 = arctan(/>i/</ 2 )- Local variations 
of the inflaton decay rate, T, produce a primordial density perturbation {Sj 


c 


1 JT 

6Y ' 


(49) 


We require T^/T 3> e^ 1 to produce a nearly scale-invariant spectrum of density perturbations that dominates over 
the inflaton perturbations. However as all the fields are oscillating on the same timescale, to - 1 , at the end of inflation, 
due to the rotational symmetry of the potential, it is not clear that one can consistently implement the modulated 
reheating proposal of Refs. j|f9] which assumes that the decay rate depends on the VEV of some light modulus field, 
which is assumed to be approximately constant during the decay of the inflaton. 

One possibility is a curvaton mechanism where the real and imaginary parts of the complex field decayed at 
different rates, Ti and IV However it turns out that perturbations from a curvaton mechanism in this model are 
always sub-dominant with respect to those from inflation perturbations. The maximum expansion from the end 
of inflation to the decay of the curvaton is N = (2/3) ln(U e /Ti) < 60, assuming the real part, say, decays before 
primordial nucleosynthesis, compared with N = (1/2) ln(7L e /Ti) to reach the same radiation density if the imaginary 
part decayed immediately at the end of inflation, T 2 ~ H e . Assuming N is a smooth monotonic function of the phase, 
0, we can estimate the perturbed expansion due to fluctuations in the phase as 


From Eq. 0 we thus have 


c 

OO Z7T <T 


^ 75 

Vciso K 4# 



(50) 


(51) 


Comparing this with Eq. 11351 1 we see that the density perturbations from the curvaton are always subdominant with 
respect to those from inflaton perturbations. 

Instead we considered a model of preheating at the end of inflation where particle production due to the resonant 
amplification of vacuum fluctuations is sensitive to the phase of a complex scalar field. In this case we showed 
that the isocurvature field fluctuations during inflation can perturb the local number density of particles produced 
in the preheating field, which leads to a primordial density perturbation. Combining the results from the inflaton 
perturbations (l35l) and the isocurvature perturbations m we obtain 





N + 


27rcr 2 


2 

k—aH 


(tan 0(1 — I? 2 )\ 2 \ 
V R 2 + tan 2 6 ) ) 


(52) 


We see that if R = gi/g 2 1, corresponding to an interaction that strongly breaks the rotational symmetry, then 
the primordial density perturbation can indeed be dominated by the nearly scale-invariant spectrum of isocurvature 
field fluctuations. Even in this model it is only possible along classical trajectories where tan 9 <C 1. 

One should bear in mind that lowering the energy scale during inflation, set by the mass to, also reduces N, the 
number of e-foldings between our present horizon scale crossing outside the Hubble scale during inflation and the end 
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of inflation 0 . Throughout this paper we have taken N ~ 60 when giving numerical values which is nearly the 
maximum value possible. Reducing N also reduces the relative contribution of inflaton perturbations in Eq. P i. 

Our model of asymmetric instant reheating is very similar in spirit to that of Kolb et al 130. Their model was 
based on a symmetric preheating interaction following inflation driven by a quadratic potential with broken symmetry. 
The spectral index of the isocurvature perturbations during inflation in that case can deviate from scale-invariance 
due to the broken symmetry. 

Ultimately, observational data will determine whether the primordial density perturbation has a significant devi¬ 
ation from scale-invariance, as predicted by adiabatic fluctuations in the inflaton field driving chaotic inflation, or 
is compatible with the very nearly scale-invariant spectrum that could be produced from isocurvature perturbations 
during chaotic inflation. An equally important ingredient in the observational constraints is the amplitude of the pri¬ 
mordial gravitational wave background that is predicted. The amplitude of gravitational waves is determined directly 
by the energy scale during inflation and is given by m 


32 

37r \M P J 

Adopting the WMAP normalisation for Vq ~ 2 x 10 -9 j2(|, the primordial tensor-scalar ratio is given by 


N - 


27TCTg 

«I7 


(53) 


'Pt i 
r = —^- ~ lO" 1 

n 


10 _6 M, 


PI 


(54) 


Hence we see that the tensor-scalar ratio is reduced if m < 10 ~ 6 Mpi. Thus if isocurvature field fluctuations contribute 
to the primordial density perturbation t he g ravitational wave background is smaller than if the density perturbations 
are due solely to inflaton perturbations (nj, Qjl O ■ 

If observations rule out the relatively large tensor contribution, and deviation from scale-invariance of the scalar per¬ 
turbations, then single-field chaotic inflation models may be ruled out. But if the primordial perturbations originated 
from isocurvature field perturbations during chaotic inflation driven by a symmetric potential, then the primordial 
perturbations could be very close to a Harrison-Zel’dovich spectrum, with n = 0.9998, and a much smaller amplitude 
of gravitational waves. 
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APPENDIX A: HANKEL FUNCTION METHOD TO CALCULATE THE ISOCURVATURE TILT 


The equation of motion for the isocurvature perturbation is 

Sx + 3 HS\ + + m 2 ^J 5x = 0. 

This can be rewritten in conformal time and using v = aSx as 

k 2 + m 2 a 2 — v = 0. 



(Al) 


(A2) 


To write the equation above in a more useful form two useful equations, that can be derived using the slow-roll 
approximation, are r = — (l/aH)(l + e + 7/3e 2 + 0(e 3 )) and V/3M 2 H 2 = 1 — en/3 where e# = a 2 /2M 2 H 2 = 
e — 2/3e 2 + 0(e 3 ). It then follows that 

m 2 a 2 = — e(3 + 5e + 0(e 2 )) , (A3) 

— = 4r(2 + 3e + 10 e 2 + 0(e 3 ), (A4) 

a t z 

and hence 


v" + 




3 6 


0{e 


v = 0 where 


v = 


(A5) 
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This equation can be solved in terms of Hankel functions assuming v is constant. Note that within a few Hubble 
times v is constant up to third order in slow roll, the solution is thus 


d(k)Hl 1) (~kr) + c 2 (k)Hl 2 '> (-kr) 


(A6) 


Matching this to the small scale Klein Gordon solution v = e lkT /V2k implies c\(k) = e *("+i/2)7r/2^jy2 and c 2 (k) = 0. 
The behavior of the Hankel function for small argument, z, is 





l/ty-K/‘2i2V—3/2 


n-) 

r(3/2) 


(A7) 


Hence the late time solution is 

|H^^(-fcr)^(l + 0(e 2 )) (A8) 

and the power spectrum is 




27t 2 \ a 


1 (-r) 1 11 


(27T) 2 


-k 3 ~ 2l/ (l + 0(e 2 )). 


(A9) 


From this form of the power spectrum it is simple to read off the spectral index to second-order in slow-roll parameters: 

10 


An v = 3 — 2 1 / ~- e 1 . 


(A10) 
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